Scaling relations for magnetic nanoparticles 
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A detailed investigation of the scaling relations recently proposed by d'Albuquerque e Castro et 
al^ to study the magnetic properties of nanoparticles is presented. Analytical expressions for the 
total energy of three characteristic internal configurations of the particles are obtained, in terms of 
which the behavior of the magnetic phase diagram for those particles upon scaling of the exchange 
interaction is discussed. The exponent r\ in scaling relations is shown to be dependent on the 
geometry of the vortex core, and results for specific cases are presented. 



I. INTRODUCTION 

In recent years, a great deal of attention has been fo- 
cused on the study of regular arrays of magnetic particles 
produced by nano-imprint lithography. Besides the ba- 
sic scientific interest in the magnetic properties of these 
systems, there is evidence that they might be used in the 
production of new magnetic devices, or as media for high 
density magnetic recording 2 . One of the main points in 
the study of such systems concerns the internal magnetic 
structure of the nanoparticles as a function of their shape 
and size. For example, in the case of cylindrically shaped 
particles produced by electrodeposition, the internal ar- 
rangements of the magnetic moments have been identi- 
fied as been close to one of the following three (idealized) 
characteristic configurations, namely ferromagnetic with 
the magnetization parallel to the basis of the cylinder 
(Fl), ferromagnetic with the magnetization parallel to 
the cylinder axis (F2), and a vortex state, in which most 
of the magnetic moments lie parallel to the basis of the 
cylinder (Vf®. The occurrence of each of these con- 
figurations depends on geometrical factors, such as the 
linear dimensions of the cylinders and their aspect ratio. 
Clearly, for the development of magnetic devices based on 
those arrays, knowledge of the internal magnetic struc- 
ture of the particles is of fundamental importance. 

Experimentally, attempts have been made to deter- 
mine, from the analysis of hysteresis curves^', the range 
of values of diameter D and height H of cylindrically 
shaped particles for which the internal arrangement of 
the magnetic moments is close to either one of the two 
ferromagnetic configurations (Fl or F2) or to the vortex 
one (V). However, such approach does not allow a clear 
description of the magnetic structure of individual cylin- 
ders, since in many cases the internal magnetic config- 
urations are not readily identifiable from magnetization 
curves. 



On the other hand, theoretical determination of the 
configuration of lowest energy of particles in the size 
range of those currently produced, based on a micro- 
scopic approach and using present standard computa- 
tional facilities, is out of reach. The reason is the exceed- 
ingly large number of magnetic moments within such par- 
ticles, which may exceed 10 9 . Recently, d'Albuquerque e 
Castro et aZp^have proposed a scaling technique for deter- 
mining the phase diagram giving the configuration of low- 
est energy among the three above mentioned characteris- 
tic magnetic configurations. They have shown that such 
diagram can be obtained from those for much smaller 
particles, in which the exchange interaction J has been 
scaled down by a factor x < 1, i.e. for J' = xJ. The 
diagram for the full strength of the exchange interaction 
is then obtained by scaling up the D' and H' axes in the 
phase diagram for J' by a factor l/x v . In their work, the 
exponent 77 has been determined numerically from the 
position, as a function of x, of a triple point (D t ,H t ) in 
the phase diagram where the three configurations have 
equal energy. The scaling technique has been applied to 
the determination of the phase diagram of cylindrically 
shapecP and truncated conicaP particles. In both cases, 
r\ turned out to be approximately equal to 0.55. 

We recall that the vortex configuration exhibits a core 
region within which the magnetic moments have a non- 
zero component parallel to the axis of either the cylinder 
or the truncated cone. We remark that the determina- 
tion of the geometry of the core (i.e. its shape and size), 
on the basis of a microscopic model in which the individ- 
ual magnetic moments are considered, would require a 
prohibitively large computational effort. For this reason, 
d'Albuquerque e Castro et aL^and Escrig et al^ adopted 
a simplified representation of the vortex core, consisting 
of a single line of magnetic moments along the axis of 
either the cylinders or the truncated cones. The phase 
diagrams thus obtained are in good agreement with ex- 
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perimental data, provided appropriate values of the ex- 
change are considered. 

The scaling technique represents a useful tool for 
studying the magnetic properties of nanosized particles. 
It is conceptually simple and rather interesting from the 
theoretical point of view. Its implementation depends on 
the determination of the exponent 77 in the scaling factor, 
which so far has been done numerically. The agreement, 
within error bars, between the values of 77 for cylinders 
and truncated conical particles suggests that this param- 
eter does not depend on the shape of the particles. How- 
ever, there still remains the question regarding the possi- 
ble dependence of 7/ on the geometry of the vortex core. 
The present work aims precisely at clarifying this point. 

We focus on cylindrically shaped particles, for which a 
large amount of experimental data is available. We adopt 
a continuous model for the internal magnetic structure of 
the particles, on the basis of which analytical results for 
the total energy in each configuration can be obtained. 
We use these results to investigate the behavior of the 
phase diagrams under scaling transformation, from which 
the value of 77 can be determined. We find that the value 
of 77 does depend on the geometry of the vortex core. 
This point is discussed at length below. 



II. CONTINUOUS MAGNETIZATION MODEL 

We adopt a simplified description of the system, in 
which the discrete distribution of magnetic moments is 
replaced with a continuous one, defined by a function 
M{r) such that M(r) SV gives the total magnetization 
within the element of volume SV centered at r. This 
model provides a fairly good basis for the discussion 
of the magnetic properties of nanosized particles. For 
cylindrically shaped particles, the magnetization density 
M(r) in the two ferromagnetic configurations, Fl and 
F2, is given by Mqx and Moz, respectively. Here Mo is 
the saturation magnetization density, and x and z are 
unit vectors parallel to the basis and to the axis of the 
cylinders, respectively. For the vortex configuration, wc 
assume that the magnetization density has the form 



M z {p) z + M v (p) (p , 



(1) 



where I and tp are unit vectors in cylindrical coordinates, 
and M z and M v satisfy the relation M 2 + M 2 = Af 2 . 
Thus, the profile of the vortex core is fully specified by 
just giving the function M z {p). It is worth pointing out 
that the functional form in Eq.([T]) does not take into 
account the possibility of a dependence of the core shape 
on coordinate z. 

We then look at the total energy of the three config- 
urations under consideration, from which the magnetic 
phase diagram can be obtained and its behavior under 
scaling investigated. We restrict our discussion to arrays 
in which the separation between cylinders is sufficie ntly 
large for the interaction between them to be ignoredP^ 



The internal energy per unit of volume, E to t , of a single 
cylinder is given by the sum of three terms corresponding 
to the magnetostatic {E^p), the exchange (E ex ), and the 
anisotropy (Ek) contributions. However, in the case of 
particles produced by electrodeposition, the crystalline 
anisotropy term is much smaller than the other tweff^, 
so its inclusion has little effect on the phase diagram. In 
view of that, it will be neglected in our calculations. 



A. Ferromagnetic configurations 

Since the exchange term depends only on the relative 
orientation of the magnetic moments, it has the same 

(F) 

value Ee X in the two ferromagnetic configurations. Since 
it also appears as an additive term in the expression for 
exchange energy in the vortex configuration, it can be 
simply left out in our calculations. 

The magnetostatic term is generally given b}^ 
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where U ( r ) is the magnetostatic potential. In the above 
expression, an additive term independent of the configu- 
ration has been left out. For the ferromagnetic configu- 
rations we find that 
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where a — Fl, F2, and N a are the demagnetizing fac- 
tors, given in SI unities bjff^ 
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In the above two equations, 2F\{a,b 1 c,d) is a hyperge- 
ometric function. Notice that demagnetizing factors de- 
pend on just the ratio D/H. 



B. Vortex configuration 

Assuming that M(r) varies slowly on the scale of the 
lattice parameter, the exchange term for this configura- 
tion can be approximated bjffQ 



ex v 



J |(vm I ) 2 + pmjV (Vm 2 ) 2 j dV , 



where A is the exchange stiffness constant, and = 
Mi/Mo, for i = x, y, z. We recall that A is proportional 
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to the exchange interaction energy J between the mag- 
netic moments.^ Making use of the expression for M(T^) 
m Eq.Q, we find 



E^ = 



2A 



f{p) P dp 



(6) 



where R = D/2, and f(p) = (dm z /dp) 2 /(l - m*) + 
(1 - m 2 z )/p 2 , with m z (p) = M z (p)/M . The additive 

term £e X on the r.h.s. of the above equation has been 
omitted. 

The magnetostatic term can be also written in terms 
of M{r). In the vortex configuration, the magnetostatic 
potential is given by 



1 



\r - n| 4tt Js 2 |r-r 2 | 



where Si and S2 are the surfaces of the top and bottom 
basis of the cylinder, respectively. After some manipula- 
tions, the expression for U (f) reduces to 



U(p,z) 



p'dp'M z (p') 



(9) 



where 7n = \H[2n]-nH [-JL]. Here, H [k] = £* =1 V* 
are the harmonic numbers. For the dipolar energy term 
111 Eq.0 we obtain 



£(V)_6W>|, 



where 
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F(n,^)= 3 F 2 
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(13) 
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Here, 3-F2 denotes the generalized hypergeometric func- 
tion. 



dk J Q {kp) J {kp')(- er kz + e -*( ff -*)), 



IV. RESULTS 



where Jo (a?) is the cylindrical Bessel function of order 
zero. Taking this result into Eq.(|2]), we fincp2 



E S = 7 Y L dk \ IP Mkp) M,(p) dp) (1 - e~ Hk ) 



(7) 



III. TOTAL ENERGY CALCULATION AND 
SCALING TRANSFORMATION 

At this point, it is necessary to specify the function 
M z (p). Since no rigorous result regarding the shape of 
the vortex core is available, we resort to a simple but 
physically plausible approximation, given by 



M z (p) = 



_ / Mo (1 - (p/ Pc ) 2 Y 







for p < p c 
otherwise 



(8) 



where p c < R and n is a non-negative constant. Alter- 
native expressions for M z {p) have been proposed in the 
literature.^ 

The above functional form for M z (p) allows us to eval- 
uate the energy integrals in Eqs. ^ and ^ analytically. 
Then, for integer values of n, the expression for E\" 
Eq.Q reduces to 



(V) 

ex m 



Having evaluated all relevant contributions to the to- 
tal energy in the three cases of interest, we are in a posi- 
tion to investigate the magnetic phase diagram for cylin- 
ders. In particular, we can look at the position of the 
triple point (D t ,H t ) as a function of the factor x which 
scales the stiffness constant A (or exchange interaction 
J). We notice that since the energy of the two ferromag- 
netic configurations, -E^ and E^ , are equal at the 
triple point, we immediately get the equation 

whose solution is £t = D t /H t — 1.10317... (independent 
of A or J). 1 ^ As a consequence, D t and H t are propor- 
tional and must exhibit the same functional dependence 
on x (or equivalently, on A). 

We proceed in our analysis by looking at the case con- 
sidered by d'Albuquerque e Castro et al.p^ in which the 
core radius is independent of x, and of the order of the 
lattice spacing (first core model). This corresponds to 
taking the limit p c « ii ( = D t /2 in the expressions for 
the total energy. In this limit, In (R/p c ) becomes much 
larger in modulus than j n , so that the latter can be safely 
neglected in Eq.Q. Then, Eqs. (jS), and ([To]) give 
the following equation for R t 



1 , , 2 2A , R t 
-N a p M = -52 In— - 

Z JX t p c 



(14) 
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where a is either Fl or F2. Now, if we scale down the 
exchange interaction by a factor x < 1, that is to say, if 
we consider a reduced exchange stiffness A' = xA, and 
assume that R t and the new radius at the triple point R' t 
are related according to R' t = x v R t , we find 



0.56- 



2A R t 
R? 



ol _ 2v 2A ^ 



Pc tt t p c 

This expression gives us the following equation for 77 

V 



1 Rt 

p c X 2r >~ 1 - 1 



lnx . 



(15) 



It is clear from this equation that 77 must in all cases 
be greater than 0.5. It approaches this lower bound only 
when R t is much larger than the lattice spacing (i.e. 
Rt Pc), in other words, when the particles have macro- 
scopic sizes. 



The behavior of r\ in Eq.(15l is presented in Fig.(l 



Fig. (la) shows 77 as a function of Rt, for 20 nm < R t < 
100 nm, and p c = 0.2 nm. We notice that in this range 
of Rt, 0.54 < 77 < 0.58. It is also interesting to look at 
the behavior of 77 as a function of x. Fig. (lb) shows 77 
as a function of x, for 0.01 < x < 1, Rt = 44 nm, and 
p c = 0.2 nm. From the curves in Figs, (la) and (lb), we 
find that for x > 0.05, 77 turns out to be close to 0.55, as 
numerically obtained by d'Albuquerque e Castro et aiP 
It is worth commenting on the effect of using a single 
value of 77, say 0.55, to scale phase diagrams for the core 
model considered just above. As already pointed out, the 
diagram for the full strength of the exchange interaction 
can be obtained from the one corresponding to a reduced 
interaction J' = xJ (with x < 1) by multiplying the axes 
H' and D' of the latter by l/x 11 . Thus, an inaccuracy 
5i] in the value of 77 results in inaccuracies 5H and 5D 
in the coordinates in the scaled diagram. Indeed, if we 
write 77 = 770 ± 5rj, with Sr]/r]o <C 1, we immediately get 
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FIG. 1. Exponent rj given by Eq.( 15 1 plotted as a function of 



Rt, for x—0.1 (solid line) and 0.2 (dashed line) (a), and as a 
function of x, for i? t =44 nm (b). 
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where H = x^H'. Since 770 0.55 and 6i]/r) Q rj 0.01 
(estimated from Fig. (la)), we find that, even for x as 
small as 0.05, the relative error 8H/Hq is smaller than 2 
%. Thus, we do not expect large discrepancies between 
the calculated phase diagram and the experimental data 
resulting from such inaccuracy in 77 since a relative error 
of 2% should not exceed the experimental error. 

We remark that the above results for 77 hold also when 
the core radius corresponds to several interatomic dis- 
tances and is kept fixed as the exchange interaction is 
scaled up or down. 

We next consider the case in which p c is adjusted 
so as to minimize the energy of the vortex con- 
figuration (second core model). From Eqs.([9| and 



( 10 1 we obtain the following equation for p c 

3a n ^-(3 n ^F(n,p c /H) 



2(ti + 2) H 6 



G(n,p c /H) 



2.4 



p M$H 2 ■ 



where 

G(n, Pc /H) 
3F2 



|, 2,|+n|, {3 + 77,4 + 277},-^! 



Eq. (??) can be solved numerically for p c in terms 
of H, A, and 77. We remark that for the core model 
under consideration, p c does not depend on the radius 
R. This follows from the fact that the outer region of 
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the cylinder does not interact with the core (apart from 
the exchange interaction across the interface between the 
two regions). As a consequence, for a given value of p c , 
the difference between the total energy of two cylinders 
of the same height but different radii docs not depend 
on p c , hence it does not contribute to the derivative of 
with respect to p c . That is to say, the equation 
for p c which minimizes the total energy of the vortex 
configuration is independent of R. 

Figure (2) illustrates M z (p) for A = 87.39 meV/nm, 
Mq = 1.4 x 10 6 A/m, and two values of H, namely 20 
and 100 nm. For each H, results are presented for n = 
2 (dotted line), 4 (dashed line), and 10 (solid line). The 
values of A and Mo correspond to those for Co, and have 
been taken from Ref.[14]. The value of p c in each case 
has been obtained from Eq.(??). 




10 p[nm] 15 



FIG. 2. Reduced magnetization m z = M z /Mo as a function 
of p, for n — 2 (dotted line), 4 (dashed line), and 10 (solid 
line). The two sets of curves correspond to H = 20 nm, and 
100 nm. Values of A and Mq have been taken from Ref. [14], 
and correspond to those for Co. 

In order to investigate the behavior of magnetic phase 
diagram upon scaling of the exchange interaction for this 
second core model, we take n = 4, which according to 
Fig. (2) provides a physically sound description of the 
core profile, and calculate the phase diagrams for dis- 
tinct values of x. Fig. (3) shows results for cylinders of 
Co (A = 87.39 meV/nm and M = 1.4 x 10 6 A/m) corre- 
sponding to x — 0.12 (dashed lines) and x = 24 (dotted 
lines). 

We then find that, for the present core model, the co- 
ordinates (D t , H t ) of the triple point follow the relations 

D t {x) = 25.61 a; 05 (16) 
H t {x) = 23.22 x - 5 (17) 

in which rj = 0.5. The diagram for the the full strength 
of the exchange interaction, x = 1, is represented by solid 
line. 

We remark that these results holds for any other inte- 
ger values of n, the reason being the fact that since p c 
is adjusted so as to minimize the energy in the vortex 
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FIG. 3. Phase diagram for Co cylinders corresponding to 
x = 0.12 (dashed lines), x = 0.24 (dotted lines) and x = 1 
(solid lines) obtained for core model 2 (see text). 



configuration, the effective radius of the core turns out 
to be independent of n, as clearly shown in Fig. (2). 



V. CONCLUSIONS 

We have carried out a detailed analysis of scaling tech- 
nique recently proposed by d'Albuquerque e Castro et at 1 
to investigate the magnetic phase diagram of nanoparti- 
cles. As already pointed out, this technique enables us 
to obtain the phase diagram for particles in the nanome- 
ter size range from those corresponding to much smaller 
particles, in which the exchange interaction has been re- 
duced. The scaling technique is easily implemented and 
represents a rather useful tool for dealing with nanoparti- 
cle systems. In addition, the existence of scaling relations 
and their connection with the model adopted to describe 
the magnetic particles bring about interesting theoretical 
considerations. 

The present work sheds light on a very interesting fea- 
ture of the scaling relations, namely the dependence of 
the exponent 77 on the model adopted for describing the 
core of the vortex configuration. Based on a continu- 
ous magnetization model, we were able to derive analyti- 
cal expressions for the total energy in each configuration, 
which allowed us to determine the exponent rj. We found 
that in the case of nanoparticles for which the core di- 
mensions, and consequently its contribution to the total 
energy, can be either neglected or do not change much 
upon scaling of A, rj turns out to be weakly dependent 
on x and quite close to 0.55. Nevertheless, when the 
contribution from the core is relevant and its size upon 
scaling of A changes so as to minimize the total energy in 
the vortex configuration, then rj becomes exactly equal 
to 0.5. 
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